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_j Abstract 

^4— > 
C/3 ' We study the dynamics of a tracer particle, which performs a totaUy directed random walk 

in an adsorbed monolayer composed of mobile hard-core particles undergoing continuous ex- 
changes with a vapour phase. In terms of a mean-field-type approach, based on the decoupling 
of the tracer-particle-particle correlation functions into the product of pairwise, tracer-particle 
' O ' correlations, we determine the density profiles of the monolayer particles, as seen from the 

i^ , stationary moving tracer, and calculate its terminal velocity, Vtr- In the general case the latter 

fi ' is determined implicitly, as the solution of a certain transcendental equation. In two extreme 

limits of slow and fast monolayer particles diff'usion, we obtain explicit asymptotic forms of 
Vtr- We show next that the density profile in the monolayer is strongly inhomogeneous: In 
front of the stationary moving tracer the local density is higher than the average value, pL, 
" ' . and approaches pL as an exponential function of the distance from the tracer. Past the tracer 

^1^ ' the local density is lower than p^ and the approach to p^ may proceed differently depending 

^O ' whether the particles number in the monolayer is not or is explicitly conserved. In the former 

■^ , case the approach is described by an exponential dependence with a different characteristic 

r^v ' length, compared to the behavior in front of the tracer; in the latter case, the density tends 

to PL algebraically. The characteristic lengths and the amplitudes of the density relaxation 
functions are also determined explicitly. 



'"^ ' Key words: Hard-core lattice gas, Langmuir adsorption/desorption model, biased tracer diffusion. 
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1 Introduction 

When a gas or vapour is brought into contact with a clean sohd surface, parts of its molecules 
will become reversibly attached to the surface in the form of an adsorbed layer. Knowledge 
of equilibrium and dynamical properties of such adsorbed layers is essential for many practical 
applications, including coating, gluing and lubrication. 

Following the seminal work of Langmuir (see, e.g. Ref.|l|]), the equilibrium properties of 
adsorbed layers have been extensively studied and a great number of important developments 
have been made. In particular, further analysis included more realistic forms of intermolecular 
interactions or allowed for the possibility of multilayer formation. As a result, different possible 
phase transformations have been predicted and different forms of adsorption isotherms have been 
found, which well explain available experimental data (see, e.g. Refs.[|]-^). 

Some effort has been also invested to the understanding of molecular diffusion in adsorbed 
layers, which has a strong impact on their global dynamical behavior. For instance, diffusional 
processes control the rates of spreading of molecular films on solid surfaces |Q,|^], spontaneous or 
forced dewetting of monolayers &fl] or island formation [ pO| . Here, some approximate analytical 
results have been obtained for both dynamics of an isolated adatom on a corrugated surface and 
collective diffusion, describing spreading of the macroscopic density fluctuations in interacting 
adsorbates being in contact with the vapour phase O-n. On the other hand, most of available 
studies of the tracer diffusion in adsorbed layers, which provides a useful information about their 
intrinsic viscosity, pertain to strictly two-dimensional models excluding the possibility of particles 
adsorption and desorption (see, e.g. Refs.[p^-pO[| and references therein). Except for Ref.Q, 
which studied driven tracer dynamics in a somewhat artificial one-dimensional model, analysis of 
the tracer diffusion in adsorbed monolayers undergoing exchanges with the vapour phase seems to 
be lacking at present. 

In this paper we study the diffusion properties of a driven tracer in a two-dimensional adsorbed 
monolayer undergoing continuous exchanges with the vapour. The system we consider consists of 
(i) a solid substrate, which is modeled in a usual fashion as a regular square lattice of adsorption 
sites supporting at most a single occupancy, (ii) a monolayer of adsorbed, mobile hard-core parti- 
cles in contact with a vapour and (iii) a single, hard-core tracer particle. Each monolayer particle 
moves randomly along the lattice by performing hopping motion on the vacant neighboring lattice 
sites. The monolayer particles undergo continuous exchanges with a vapour phase, i.e. may desorb 
from and adsorb onto the lattice with some prescribed rates dependent on the vapour pressure, 
temperature and the interactions with the solid substrate. On the contrary, the tracer particle 
is constrained to move on the two-dimensional lattice only, i.e. it can not desorb to the vapour, 
and it is subject to a constant external force E. Hence, the tracer performs a biased random 
walk, which is constrained by the hard-core interactions with the monolayer particles, and always 
remains within the monolayer, probing its response on the internal perturbancy or, in other words, 
its frictional properties. Here we focus on the limit of sufficiently large external force, such that 
the tracer may move only in one direction, but not large enough to make it slide regardless of the 
surface corrugation. The results for the general case of arbitrary E, which allow us to deduce the 
analog of the Stokes formula for 2D adsorbed monolayers and to define the corresponding friction 
coefficient, will be presented elsewhere p^ . 

In terms of a me an- field- type approximation of Ref. |23|| , which is based on the decoupling 
of the tracer-particle-particle correlation functions into the product of pair- wise correlations, we 



define the stationary density profiles of the monolayer particles, as seen from the moving tracer, 
and we determine analytically the terminal velocity Vtr of the tracer. We show that the adsorbed 
monolayer particles tend to accumulate in front of the driven tracer, creating a sort of a "traffic 
jam" , which impedes its motion. Thus the density profile around the tracer is highly inhomoge- 
neous: the local density of the monolayer particles in front of the tracer is higher than the average 
and approaches the average value as an exponential function of the distance from the tracer. The 
characteristic length and the amplitude of the density relaxation function are calculated explicitly. 
On the other hand, past the tracer the local density is lower than the average; we show that de- 
pending on the condition whether the number of particles in the monolayer is explicitly conserved 
or not, the local density past the tracer may tend to the average value either as an exponential 
or as an algebraic function of the distance, revealing in the latter case especially strong memory 
effects and strong correlations between the particle distribution in the monolayer and the tracer 
position. 

Further on, we find that the terminal velocity of the tracer particle depends explicitly on the 
excess density in the "jammed" region in front of the tracer. This excess density, in turn, depends 
on the magnitude of the velocity, as well as on the rate of the adsorption/desorption processes 
and on the rate at which the particles can diffuse away of the tracer. The interplay between the 
jamming of the monolayer, produced by the tracer particle, and the rate of its homogenization due 
to diffusion and adsorption/desorption processes, manifests itself as a medium-induced frictional 
force exerted on the tracer, whose magnitude depends on the tracer velocity. As a consequence 
of such a non-linear coupling, in the general case, (i.e. for arbitrary adsorption/desorption rates 
and particles diffusion coefficient), Vtr can be found only implicitly, as the solution of a certain 
non-linear equation relating Vtr to the system parameters. This equation simplifies considerably in 
the limit of small or large particles diffusivity, in which two cases exlicit asymptotic expressions for 
the tracer velocity are obtained. We finally remark, that a qualitatively similar physical effect was 
predicted recently for a different model system involving a charged particle moving at a constant 
speed at a small distance above the surface of an incompressible, infinitely deep liquid. It has been 
shown in Refs.[g4[g^, that the interactions between the moving particle and the fiuid molecules 
induce an effective frictional force exerted on the particle, producing a local distortion of the liquid 
interface, - a bump, which travels together with the particle and increases effectively its mass. The 
mass of the bump, which is analogous to the jammed region appearing in our model, depends itself 
on the particle's velocity resulting in a non-linear coupling between the medium-induced frictional 
force exerted on the particle and its velocity |24, 25|. 



The paper is structured as follows: In Section 2, we formulate the model and introduce basic 
notations. In Section 3, we write down the dynamical equations which govern the time evolution 
of the monolayer particles and of the tracer. Section 4 is devoted to the analytical solution of 
these evolution equations in the limit i — > oo ; here we also present some general results on the 
form of the density profiles around stationary moving tracer and on the tracer terminal velocity, 
which is given implicitly, as the solution of a transcendental equation. In Section 5, we derive 
explicit asymptotic results for the tracer terminal velocity in the limits of small and large particles 
diffusivities. As well, here we discuss the forms of the density profiles at intermediate scales. 
Corresponding asymptotical behavior of the density profiles at large distances from the tracer is 
discussed in Section 6. Finally, we conclude in Section 7 with a brief summary and discussion of 
our results. 



2 The model 

Consider a two-dimensional square lattice of adsorption sites {X, Y) , which is brought in contact 
with a reservoir of identic, electrically neutral particles (vapour phase) (Fig.l), maintained at a 
constant pressure. We suppose that the particles may leave the reservoir and adsorb onto any 
vacant lattice site at a fixed rate //t*. Then, the adsorbed particles may move randomly on 
the lattice by hopping with a rate I/At* to any of four neighboring lattice sites, which process 
is constrained by hard-core exclusion preventing multiple occupancy of any site. They can also 
desorb from the lattice back to the reservoir at a site- and environment-independent rate g/r* . We 
assume, for simplicity of exposition, that typical adsorption, desorption and jump times are equal 
to each otherR and are denoted by r*. Consequently, at any time instant the variable rj{X^Y), 
characterizing the occupation of the site {X, Y) such that 



77(x,r) 



1, if the site (X, Y) is occupied by an adsorbed particle 
0, if the site (X, Y) is empty. 



can change its value due to adsorption, desorption and constrained random hopping events. Note 
that the total number of particles in the adsorbed monolayer is not conserved in this dynamical 
process. However, the mean density pl of the adsorbate, p^ —< vi^j Y) >, approaches as i — + oo 
a constant value 

which equation is usually referred to as the Langmuir adsorption isotherm |l|| . 

Further on, at i = we introduce at the lattice origin an extra hard-core particle, whose motion 
we would like to follow and which will be called in the remainder as the tracer. We stipulate that 
the tracer is different from the adsorbed particles in that it can not desorb from the lattice and in 
that it is subject to some external driving force, which favors its jumps to a preferential direction. 
Such a situation may be realized, for instance, if this only particle is charged and the system is 
subject to a uniform electric field. Here we will focus on the particular case when the external field, 
oriented in the positive X-direction, is sufficiently strong, such that the tracer performs totally 
directed random walk (Poisson process) along the X-axis in the two-dimensional monolayer of 
adsorbed mobile, hard-core particles. The tracer can be thought off as a certain external probe, 
which is designed to measure the resistance offered by the monolayer particles to the external 
perturbance. 

More precisely, we define the tracer dynamics as follows: We suppose that the tracer, which 
is at the site Xtr{t) at time t, waits an exponentially distributed time with meann r, and then 
attempts to hop onto the neighboring site Xtj.{t) + a, a being the lattice spacing. The jump is 
actually fulfilled only if the target site is vacant at this moment of time; otherwise, i.e., if the 

^These times, T^di ''"des *nd Tjumpi respectively, can be readily made different and restored in our final results 
by the mere replacement / — > t* f/T^j_, g -^ t^q/t^^^ and I -^ It* /Tj^mp- 

■^We suppose that in the general case this mean time r is different from the corresponding time r* associated 
with the monolayer particles dynamics. As a matter of fact, this should be the case merely because the tracer- 
substrate interactions may be different from the particle-substrate ones. Varying r we can mimic different possible 
situations; in particular, t = corresponds to the case when the tracer simply slides along the substrate regardless 
of the surface corrugation. 



target site is occupied by any adsorbed particle, the jump is rejected and the tracer remains at its 
position. 

3 Evolution equations 

Let P{Xtr,ri;t) denote the probabihty of finding at time t the tracer at the site {Xtr,0) and all 
adsorbed particles in the configuration 77 = {r]{X,Y)}. Further more, let rj^-^'^^'^-^ '^ ^ denote the 
configuration obtained from 77 by exchanging the occupation variables of the two neighboring sites 
{X,Y) and {X',Y'), which describes the Kawasaki-type particle- vacancy exchange between the 
sites {X,Y) and {X',Y'). Next, we denote as r]^-^'^^ the configuration obtained from the original 
configuration 77 by the replacement ri{X, Y) ^ 1 — ri{X, Y) , which corresponds to the Glauber-type 
flip of the occupation variable due to the adsorption/desorption events. 

Then, counting up all events which can result in the configuration {Xtr, rf) at time t or modify it, 
we write down the following master equation, which governs the time evolution of the configuration 
probability P{Xtr, tj; t): 

PiXtr,mt) = ^ ^ "{P(X„,,ry(^.^)^(^+-.^);t)-P(X„,,r;;t)} + 

^ {X,Y) 



+ -{(1 - r^Xtr. 0))P{Xtr - a, 77; t)-{l- r^iXtr + o, 0))P(Xt,, 77; i)} + 

r 

+ 4 

T* 

{X,Y) 



J2'{{1- v{X, Yj)P{X,r, 77(^-^); t) - V{X, Y)P{Xtr, m t)} 



4 E 'MX,Y)P{Xtr,v^'''''^;t)-{l-ijiX,Y))P{Xtr,mt)h 

^ (X,Y) 



(2) 



where the dot denotes the time derivative, double prime indicates that the summation extends 
over all lattice sites {X,Y), excluding the special sites {Xtr — tr, 0) and {Xtr,0); a single prime 
after the summation sign means that the sums runs over all lattice sites excluding the site {Xtr , 0) 
only. 

Now, we can obtain the instantaneous velocity Vtr (t) of the tracer particle directly from Eq. (0) . 
To do this, we multiply both sides of Eq.(||) by Xtr and sum over all possible configurations {Xtr, ??), 
which yields the following equation for the instantaneous tracer velocity 



Vtr{t)^^^ = ^{l-k{aA,t)), (3) 

where the bar denotes the configurational average and 

fc(Ai,A2;i)= Y. ri{Xtr + \i,\2)P{Xtr,mt) (4) 

{Xtr,r)) 

is the probability of having at time t an adsorbed particle at position (Ai, A2), defined in the frame 
of reference moving with the tracer particle. In other words, A;(Ai, A2) can be interpreted as the 
density profile as seen from the moving tracer. 



Consequently, Eq.p) relates the instantaneous velocity of the tracer particle to the particle 
density in the immediate vicinity of the tracer. Note, that if the monolayer was perfectly stirred, 
i.e. fc(Ai, A2) = PL everywhere, (which implies decoupling of Xtr and 77), we should obtain from 
Eq.(0) a trivial mean- field result 

V,r^^-^^^J^, (5) 

T 

which states that the tracer jump time t is merely renormalized by a factor (1 — p^)^^, which is 
the inverse concentration of voids in the monolayer: (1 — Pl)/t defines simply the frequency of 
successful jump events. However, this turns out not to be the case and fc(Ai, A2) appears to be 
different of the equilibrium value pL everywhere, except for infinitely separated sites, A1.2 — » ±cx3. 
Moreover, as we proceed to show, fc(Ai, A2) is itself dependent on the tracer velocity which yields 
ultimately a strongly non- linear equation determining Vtr- 

Hence, in order to calculate the mean velocity of the tracer we have to determine the form of the 
density profile fc(Ai, A2), or more precisely, the value of the mean density at the neighboring to the 
tracer site (Ai = cr, A2 = 0). The latter can be found from the master equation (0) by multiplying 
both sides of Eq.(0) by r]{Xtr + Ai; A2) and performing the summation over all configurations. 

In doing so, we find the following set of equations, which hold for the sites separated from the 
tracer by the distance exceeding the lattice spacing a, i.e. such that A^ + Aj > (J^, 

+ - y {l-f]{Xtr+a,0))P{Xtr,mt){v{Xtr + M+T,X2)-v{Xtr+Xi,X2)}, (6) 

where the symbols A^j 2 denote the central second-order finite difference operators of step a, 

AAifc(Ai, A2; t) = fc(Ai + a, A2; t) + fc(Ai - cr, A2; i) - 2fc(Ai, A2; t) (7) 

On the other hand, for the sites adjacent to the tracer position, i.e. ((7,0), (— cr, 0) and (0, ±ct), 
respectively, we get 

At* t* t* 

+ - y [1 - viXtr + <7,0))P{Xtr,v;tHXtr + 2a,0), (8) 

T ^ — ' 
(Xtr-.ri) 

k{-a,0;t) = ^{Va, +AAjfc(-a,0;t)-^^fc(-a,0;t) + 4 + 
4t-* y* t* 

- y {l-l]{Xtr + CT,0))P{Xtr,V;tHXtr-a,0), (9) 



'^ {Xtr,V) 



and 



fc(0, ±a; t) = /^{Va, + Aa, }fc(0, ±a; t) - ^k{0, ±a; t) + ^ + 

- y {l-r]iXtr + c7,0))P{Xtr,v;t){v{Xtr + T,±a)-j^{Xtr,±a)}, (10) 



'^ {Xt^,v) 



where Vai 2 stand for the forward difference operators, of the form 

VA,fc(Ai,A2;t) = fc(Ai + (T,A2;i)-fc(Ai,A2;t), (11) 

and 

VA2fc(Ai,A2;t) = fc(Ai,A2 + (T;i)-fc(Ai,A2;i), (12) 

Note, that Eqs.(g) to ( [lO| ) are different from Eq.(|g), since the density profile at the neighboring 
to the tracer sites is strongly perturbed by its asymmetric hopping rules. As a matter of fact, 
Eqs.(H) to ( pi]| ) can be thought off as the boundary conditions for Eq.(g). 

Now, several remarks on Eqs.(^) to ( [lO| ) are in order. First of all, one notices that these 
equations are not closed with respect to fc(Ai,A2) but rather coupled to the third-order, tracer- 
particle-particle correlation functions. In turn, if we attempt to derive the evolution equations for 
the third-order correlation functions, we find that the latter appear to be coupled to the fourth- 
order correlations. Consequently, in order to determine the tracer velocity, one faces the problem 
of solving an infinite hierarchy of coupled equations for the correlation functions. Here we resort 
to the simplest non-trivial closure of Eqs.(M) to ( pO| ) in terms of the pairwise correlation function 
A:(Ai, A2), which is based on the following decoupling approximation of the third-order correlation 
functions. 



J2 V{Xtr+Xl,X2){l-v{Xtr + 'J,0))P{Xtr,V;t)~ 

~ { Y. ^(^tr + Ai,A2)F(Xi„77;i)}x{ ^ {1 - ri{Xtr + a,0))P{Xtr,mt)} = 

{Xt.r,ri) {Xt,.,v) 

= fc(Ai,A2;t)(l-fc(a,0;i)), (13) 

i.e. the average with the weight P{Xtr,iT,t) of the product of several occupation numbers of 
different sites is set equal to the product of their average values with the weight P{Xtr, tj; t). 

We hasten to remark that the approximate closure of the evolution equations in Eq.|l3|), as well 
as some similar but not obviously equivalent approximations [16, 17|, have been already employed 



for studying related models of tracer diffusion in hard-core lattices gases and shown to provide 
quite accurate description of the dynamical and stationary behavior. The decoupling in Eq.(|l3|) 
has been first introduced in Ref . |23t| to determine the properties of the driven tracer diffusion 
in a one-dimensional hard-core lattice gas with the conserved number of particles, i.e. without 
exchanges of particles with the reservoir. Extensive numerical simulations performed in Ref. 123] 
have demonstrated that such a decoupling is quite a plausible approximation for the model under 
study. Moreover, rigorous probabilistic analysis of Ref.pO] has shown that for this model the 
results based on the decoupling scheme in Eq.(|l^) are essentially exact. Furthermore, the same 
closure procedure has been applied recently to study spreading of a hard-core lattice gas from a 
reservoir attached to one of the lattice sites |^ , and to treat the biased tracer dynamics in a one- 
dimensional model of adsorbed monolayer in contact with a vapour phase, i.e. a one-dimensional 
version of the model to be studied here. Also in these cases an excellent agreement has been 
observed between the analytical predictions and Monte Carlo simulations data. Last but not 
least, as we set out to show elsewhere pa, in case of arbitrary force exerted on the tracer the 



decoupling in Eq.(0) reproduces the results of Refs.|16|| and |T^ for the tracer diffusion coefficient 
in two-dimensional hard-core lattice gases with a conserved particles number (the limit f,g — 0, 
f/g = const in our case); these results are known to be exact in the limits of small and large 
particle densities and provide a very accurate approximation for the tracer diffusion coefhcient in 
two-dimensional hard-core lattice gases with arbitrary particle density |15[| . We thus expect that 
it will render a plausible description of the tracer dynamics in the two-dimensional model under 
study and base our further analysis on this mean-field-type approximation. 

Using the approximation in Eq.(|l3|), we can rewrite Eq.(o) in the following closed form 



fc(Ai,A2;t) = Zfc(Ai,A2;t) + 4, (14) 

r 

where the operator L is given by 

L = J-{Aa, + Aa,} + ^{1 - fc(a,0;t)}VA, - ^ (15) 

Further more, we find 

fc(cr, 0; t) = (l k{Xi, As; t)) + jL + i{i _ k{cr, 0; t)}k{a, 0; t) + -—k{a, 0; i), (16) 

V / Ai=(T,A2=o T* T At* 

fc(-a,0;t)= (L/c(Ai,A2;i)) + ^ + -^k{-a,0;t), (17) 

V / Ai=-o-,A2=o T* At* 

and 

fc(0,±cr,0:<) = fLfc(Ai,A2;t)') + — + k(0,±cr;t) (18) 

V / Ai=0,A2 = ±CT T* 4t* 

Equations (|lj) to (|l8|) constitute a closed system of equations, which suffice the computation of 
the density profiles and tracer velocity. Note, however, that these equation are non-linear, since 
/e(Ai = (7,0) enters the prefactor before the gradient term, which makes such a computation to 
be quite a non-trivial problem. Below we consider the solution of this system of equations in the 
limit t ^ oo. 



4 Stationary solution of the evolution equations 

We turn now to the limit i ^ oo and suppose that both the tracer velocity Vtr{t) and the 
density profile around the tracer attain stationary non-zero values Vtr and fc(Ai,A2), i.e. Vtr = 
limt-,ocVtr{t) and fc(Ai, A2) = limt^ooH^i, >^2;t). 

Next, it is expedient to rewrite Eqs(|l4|) to (|l8| ) in terms of an auxiliary function hn^m'- 



hn,m = k{n(T, ma) - pL, (19) 

which determines the local deviation of the density from the average value pL- In terms of this 



auxiliary function Eqs(|l4|) to (|8D become for (n,m) 7^ (0,0), (±1,0), (0,±1) 



v 

(1 + 'P)hn+l.,m + hn-l,m + /ln,m+l + /l«,m-l - 4(1 + — + S)hn,m = 0, (20) 

while hn,m ill the immediate vicinity of the tracer obey 

(1 + r)h2.o + hi,i + /ii,_i - (3 + 45)/ii,o + VpL = 0, (21) 

7i_2,o - (1 + V)h^ifl + /i_i4 + h^i^^i - (3 + 4(5)/i_i,o - T'pL = 0, (22) 

(1 + r)hi,i + h-1,1 + ho,2 - (3 + AS)ho,i = 0, (23) 



(1 + P)h,,^i + ;i_i,_i + /i,o,-2 - (3 + 4(5);io,_i = 0, (24) 

where Dq denotes the diffusion coefficient of an isolated adsorbed particle: Dq = la'^/Ar*, the 
characteristic parameter V — Vtr(y/Do is akin to the so-called " Peclet " number in the hydrody- 
namics, and we have the parameter S = (T^{f + g)/4DQT* . The parameter S can be also written as 
i5 = (/ -I- g)/l = (1 — pLJ^^ig/l), where the factor (g/l) is an exponential of the difference of the 
energy barriers against lateral diffusion and desorption. Hence, 6 compares the relative weights of 
the adsorption/desorption and the diffusion events. Note, that the factor {g/l) is usually small for 
most of realistic experimental situations, such that 5 is small provided that p^ ^ 1. On the other 
hand, at high adsorbate densities, 5 can attain relatively large values. As well, the Peclet-type 
number V can also be small or large, depending on the physical situation. Below we will consider 
different possible limits for V and S in order to elucidate the asymptotic behavior of the tracer 
velocity. 

In a standard approach, Eqs.(po|) to (|2J) can be solved by introducing the generating function 
of the local deviation hn^rm i-e. 



oo oo 

H{z,w)^ Y, Y. V™2"^«", (25) 

n— — oo m— — oo 

Multiplying Eqs.(E3) to (pih by z" and w™ and summing over all n and m, we find then that the 
generating function H{z, w) is given explicitly by 

r P (l + P)-!-! 

H{z,w)^-K{z,w)[z+iu~'iil + -+S))+ ^ ' j , (26) 

where u = w + l/w, and 

K{z, w) = ((1 + V)hi,o + Vpl) (^ - 1) + ho.i{u - 2) + (/i_i,o - ^-Oi)^^^^ (27) 



Note, that H{z, w) depends on w only in the combination w + l/w, which signifies that hn,m is, 
as it could be expected, an even function of m. 



Now, hn,m can be found directly from Eqs.(|2q) and (27) through the standard inversion for- 
mulae p7[ |, which amounts, however, to computating rather complex integrals. Here we will use 
instead a more straightforward approach, identifying ft.„.m from the expansion of H{z, w) into the 
series in Eq.pa). To do this, we first expand H{z, w) into a series in powers of (1 + 5 + 7^/4)^^, 
which gives (if z + ii^ +w<4(l + ^+(5) which is always possible) 

H{z,w)^K{z,w)Y,(^^(l + '^+5)^ (^z+^^^ + r^ (28) 

Next, expanding (z + (1 + 'P)z~^ + uf in series in powers of w, we rewrite Eg. ([2^) as a multiple 
series of the form 



where (j^) denote the binomial coefficients. Lastly, gathering the terms with equal powers of z and 
w, we find from Eq.(E3) that hn^m is given explicitly by 



— (1 + V)hi^o\/ nFn-l,rn + h-i^^V nFn,m + hos^niFn^m — V pL{Fn+l,m — ^n-l/ 



4(1 + (5 + 7^/4) 



(30) 



where 



F^^„, = (l + V)-'''^ P{n,m\0. (31) 



with 



Vl+Vt \ ^ f t 



p{nMQ - r^'^-'^wrrph))' 



{l + 5 + V/A)J "'\2{l + 5 + V/4) 



TT ^7r V i{k^n+kym) 

j2' S 



(0 \2 I I ^'k , (32) 

in which equation In{x) denotes the modified Bessel function, while the parameter C, (which 
appears to be an important control parameter and which will be repeatedly used in what follows) 
is defined as 



''"2(1 + 5 + 7^/4)' ^-^-^^ 



and 



, , _ cos{ky) + ^/\ +Vcos{h^) 

^^""^ - 1 + Vl+P ^^^ 



Note, that P{n, m\() can be thought off, in view of the form of Eq.(p2), as the generating function 
for the probabihty of occupancy of the site (n, m) (the lattice Green function) for a particular 
type of two-dimensional Polya random walk, whose structure function is defined by Eq.(|3J) (see, 
e.g. Refs.[E8| and [E3] for more details). 

Now, two comments on the just-derived expression for /i„,m, Eg. (pO[) , are in order. First of all, 
we note that Eq.(pfl) is not closed, since it contains on the right-hand-side unknown functions; 
namely, /i±i,o a-nd /io,i. Second, all hn^m are functions of the tracer terminal velocity, which still 
remains undefined. Consequently, in the remainder we have to find /i±i,o(^r) and /io,i(Vtr) as 
functions of the microscopic parameters, (such as /, g, r and etc.), and Vtr, and then, inserting 
the obtained expression for /ii,o(Vtr) to Eq.(H), derive a closed with respect to Vtr equation. 

As a matter of fact, the deviations from p^ in the immediate vicinity of the tracer, i.e. 
h±i,QiVtr) and /io,i(Vtr), can be found directly from Eq. (pO|) . Setting n = ±1, m = and 
n = 0,m = 1, we find from Eq.(po|) that h±i,Q{Vtr) and /io,i(^tr) obey the following system of 
linear algebraic equations 



A{Vtr) X h{Vtr) = VpLF{Vtr), 

where the column-matrices h{Vtr) and F{Vtr) are defined as 



h{Vtr 



while the square matrix A{Vtr) has the form 




Fo,o — F^2,o 
Fi,i — F_i_i,y 



(35) 



(36) 



-(1 + P)(ac + VnFo,o) V„Fi,o A„Fi,o 

A{Vtr)=\ -il + r)VnF-2,o -ac -I- V„F_i,o A„F_i^o 

-(1 + 7')V„F_i,i V„Fo,i -ac -t- A,„i^o,i, 



(37) 



where a^ = 2(H--\/l + 'P)/C and V„ (A^) are the finite difference operators of unit step in Eq.(0) 
((|ll[)), acting on the variable n (m). To define the tracer terminal velocity, we merely have to 
know hifi{Vtr)- It follows from Eqs. (^ to (^7|) that /ii,o(Vtr) is given by 



hi,oiVtr) = Vpi 



det hi,o(Vt,) 
detA(Vtr) ' 



(38) 



where /ii.o(Vtr) is the matrix obtained from (pTh by replacing the first column by F{Vtr)- 

Now, we are in position to obtain a closed equation, which determines implicitly the tracer 
particle terminal velocity in terms of the parameters /, g, pL and Dq. Substituting Eqs.(|38[) and 
(1^ into Eq.(|), we find that Vtr obeys 



Vtr 



^(1-Pl) 



a'^PL det hi,o(Vtr) 



Dot det A(Vtr) 



(39) 
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where the first term is the mean field result for the tracer velocity (see Eq. (m), obtained under 
assumption that the monolayer is perfectly stirred; terms in brackets stem from non-linear coop- 
erative effects associated with formation of inhomogeneous density profiles, whose characteristic 
properties depend parametrically on the tracer velocity. Equation (p9|) , which is the main general 
analytical result of our paper, is certainly too complicated to be solvable in the general case. Below 
we will study some of its asymptotical solutions. 

5 Asymptotic behavior of the tracer stationary velocity and 
density profiles at a fixed distance from the tracer 

In this section we consider the asymptotical behavior of the tracer particle velocity and density 
profiles at a fixed distance from the tracer in the limit of slow and fast particles diffusion (n and m 
are taken arbitrary but fixed, while Dq is allowed to vary). As well, we will also specify different 
regimes corresponding to the cases with V ^ 1 and V -^ 1, which may be realized at different 
monolayer densities. We note that such an analysis turns out to be rather non-trivial, which 
happens not only because of the non-linearity of the equations involved, but chiefly due to the 
non-uniformity of the asymptotical expansions with respect to the parameter S. As a matter of 
fact, such a non-uniformity is just the consequence of the qualitatively different behavior occurring 
in the system with explicitly conserved (S = 0) and non-conserved {6 > 0) particles number. 
Consequently, the results for the system without adsorption and desorption can not be obtained 
as a limit (5 ^ of the corresponding results with arbitrary (even vanishingly small) values of i5. 
Lastly, we note that the limits Dq — > oo (or Dq — > 0) and n, to — > oo do not commute; hence, 
we consider the asymptotical behavior of the density profiles at large distances from the tracer 
separately in Section 6. 

5.1 Limit of a small particles diffusivity. 

Consider first the asymptotical behavior of Vtr in the limit of a vanishingly small diffusion coef- 
ficient Do of the monolayer particles, Dq — > 0, which situation can be thought off as a weakly 
perturbed original Langmuir model (Da ~ 0). It appears, as we have already remarked, that this 
asymptotical behavior can be different in the case when the particles number in the monolayer 
is not explicitly conserved, i.e. when both /, g > 0, and in the case with a conserved particles 
number, when both f,g~0 but their ratio is fixed, f /g = pl/{^ — Pl)- Consequently, we have to 
study these two cases separately. 

A. Non-conserved particles number. We assume first that both / and g are non-zero, such that 
(5 7^ 0, which means that adsorption and desorption do change the occupation of each site. In this 
case, it seems natural to expect that the tracer will never be blocked and will always continue 
moving at a constant velocity. In other words, we anticipate that in this case the velocity attains 
a constant non-zero value V^^' > as Dq — > and hence, that V ^ 1, 5 ^ 1, which means that 
C ^ 1. We will refer to this limit oi Dq —^ Q, 5 ^ 1 and P » 1 as the limit (a). We seek then Vtr 
in the form 

Vtr = F("' + F(i) (^) + 0({^f) . (40) 
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Further on, expanding the denominator in Eq.(p3) in powers of the parameter C, we calculate 
to the first order in powers of Dqt* /a^ the functions Fn^m involved in Eq.(B9h. These are listed 



in Appendix 1. Then, substituting the obtained expressions for Fn^m into Eq.(39), we have then 



and 



„.,.(i^|,,_^j («) 



T/(i) n ^ pL + 2 + i{f + g)T/T* 



Further on, using Eqs.(El|) and (p2), we find the following explicit results for the deviations 
from the equilibrium mean density pL in the immediate vicinity of the tracer. These obey, in the 
limit (a): 



(43) 



Pl(I-Pl) L t (pi + 2 + 3(/ + .g)r/T*)/Z?or*x I ^^/,Z?or*,2 



and 

, Pl(I-Pl)' /i^oT*\ n{l^^^\A (AK^ 

As it could be expected intuitively, h\j^ > and h_iQ < 0, which means that there is a condensed, 
"jammed" region in front of the tracer, and a depleted region past the tracer. On the other hand, 
the fact that /iq.i > and that h±i^Q under certain conditions are non-monotonous functions of 
PL appear to be non-trivial. We will discuss the characteristic properties of these two regions at 
the end of this section and in the section ||. 

Note, that the result in Eq.(pl|) can be simply obtained from the initial system of equations 
( P0| ) to (M) by assuming that the deviation hnfi from the mean density is zero for n > 2, i.e. 
the density profile in front of the density profile is an abrupt step function. Then, it follows from 
Eqs.(pl|) and ( |2^ ) that /i±i,o ~ '^Vph/i^ + 4(5) Plugging this approximate result into Eq.(||) and 
turning to the limit Dq ^ we actually recover Eq.(^l]). Correction term in Eq.([4^) is thus 
associated with the appearance of a smooth, long-range density variation with the distance from 
the tracer due to diffusion, which couples effectively the evolution of /i„^m at different lattice sites. 
Note finally, that Eq.(|4l|) reduces to the essentially mean-field in Eq.(|^), i.e. Vtr — cr(l — Pl)/t^ 
in the limit r* — > 0, i.e. in the limit when adsorption/desorption processes proceed at infinitely 
fast rate resulting in a complete homogenization of the monolayer. 

Next, using Eq. (pO|), and taking advantage of the fact that to the lowest non-trivial order in 



powers of Dq the function Fn^m in Eq.(31) is given by 
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F„ 



where 



n\ J\V^°)a + (j'^{f + g)/T* 



\n\ + \m\ 



^ni 



(46) 



Jfl — 



1, if n> 0, 

(yWcr/Z^o)'"', if n<0, 



we find that, at a finite distance from the tracer, the local deviation from the equilibrium mean 
density obey: 



P(l-P) 



P+{f + g)T/T*\n-l 



-l + |m| 



D^r* p+{f + g)T/i 



^2(/ + .g)l + (/ + .g)T/T* 



2-i+k 



for n > 0, 



(47) 



and 



P(l - P) 



l + if + g)T/T*\\n\ 



1 1 + 1 m I 



1-P 



1 + (/ + 9)r/T* 



"1-1 



Dor* p+if + g)T/T* 



<^^{f + g)i + {f + g)T/T* 



ho., 



I"^|p(l-p)^ 



Dor* P+if + g)r/T* 



'J'{f + 9)l + if + 9)r/r* 



\7n\ 



for n < 0, 



(48) 



(49) 



il + {f + 9)r/T*)ip+if + g)T/T*) 

Note, that these expressions generalize the leading order contributions of Eqs.([70|)-(^5|). It may 
be also worthy to remark that hn^m turns out to be a non-separable function with respect to the 
variables n and m. The density profiles defined by Eqs. (^) to (^9|) are depicted in Fig. 2a. 

Lastly, we comment on the range of applicability of all the results of this subsection, which are 
based on the assumption that V ^ I and (^ ^ 1. This implies, in turn, that these results are only 
valid when Dq obeys the inequality 



Do<^ 



<yHi- PL) 

2t 



Plt 



if + 9)r 



(50) 



which sets up an upper bound on Dq showing how small it should be to insure the utility the 
previous expressions. 

B. Conserved particles number. Consider next the special case when both f,g —^ 0, but their 
ratio f/g is kept fixed, f/g — Pl/{^ ~ Pl), which corresponds to the situation with a conserved 
number of particles (no adsorption/desorption), i.e. a standard model of a two-dimensional hard- 
core lattice gas. In this case the parameter S is exactly equal to zero, S — 0, F^"^ should vanish 
according to Eq. (^) , and hence, Vtr is expected to scale linearly with Dq in the lowest order with 
respect to the particle diffusion coefficient. Note, however, that in the case 6 = the inequality 
in Eq.(pG) is not satisfied for any Dq except Dq = 0, which means that, in principle, Eqs.([4l|) and 
( [42| ) do not describe the behavior of the tracer velocity properly. 

Let us now examine separately the case 5 = for which the results ( [ll|) to (|45| ) do not apply. 
Assume first that the Peclet number is large, i.e. V ^ 1, which will be checked for consistency 
afterwards. Then, in this limit Dq — > 0, ^ = and V ^ 1, which will be referred to in the 
remainder as the limit (b), we have that C ^ 1 and we can use the expressions for Fn^m (see 
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Appendix 1) with S ~ 0, provided that V^^'' is replaced by Vtr- This yields for the terminal 
velocity the following result 

2(1 -Pl) 



Vtr 



crpL 



-Do 



(51) 



It follows then from Eq.(pT|) that here the Peclet number is given approximately by T' « 2(1 — 
Pl)I Pl\ hence, Eq.(|5l|) is consistent with the underlying assumption only when pi ^ 1. In this 



limit, we find that hiQ ^ 1, /i-i,o ^ ^Pl- ^nd hgi ^ Pl/'^- Note that in this case, Hq i 

is independent of Dq to leading order, contrarily to the prediction of Eq. (Ea). 

More generally, using Eq. ( |30| ) and Eqs.(p5[)- (100), we find that at a finite distance from the 
tracer 

-l+|m|\ / Q\n—l-\-\m\ 



i^n' 



.-r )(f) ^°^">«' 



and 



/l"l+l™l\ \m\ f p\\™-\ 

hn,m ~ - 1 — I — - - for n < and m yt Q 

\\n\ / |m| - n+ 1 V2/ 

ho,m - \m\ (^-j . 
Furthermore, using the expansion 

Fn,o = 1 + ^^"^ ^^ + O(^) for <5 = and n < 0, 



we have 



hn,o ~ —p for n < 



(52) 



(53) 



(54) 



(55) 



(56) 



Again, these expressions are not separable in n and m. The density profiles defined by Eqs. (|52| ) 
to ( p6| ) are depicted in Fig. 2b. 

Let us consider now the limit (c) when Do ^ 0, (5 = 0, but T-" <C 1. In this limit we evidently 
have C ~ 1. Notice now that for C = 1, the integral in Eq . (|32|) diverges at small values of the 
wave-vector k, so that when 1 — C is small but non-zero, the integral over k is dominated by 
the small-k behavior (strictly speaking, this is a priori true only if n and m are large enough 
to insure the validity of the stationary phase method) . Following MontroU and Weiss ||2^ , who 
have studied asymptotic behavior of P{n,m\C,) in Eq.(|3^) within the context of Polya random 
walk on 2D lattices, (see also Ref.|g^), wc first represent the structure function in Eq.(p4|) as 
1 - A(k) - i (alkl + alkl\ , where 



V V 



{l + Vf/^ 

(i + VTTp)i/2 



1 



{1 + ViTvy/^ 



(57) 



Then, using the stationary phase method, one finds [g9| the following representation of P[n^m\C,) 
in Eq.(|3l), which holds in the limit C ^ 1- 



P(7i,m|C) 



-A'n 



^ x^y 



1/2n 



{[^i^-oih-Ml 



(58) 
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where Ko{z) is the modified Bessel function of the third kind. In the hniit C ^ li Eq. (psh yields: 

P(n, m|C) = --ln{^) - ^"("'+"^') - ?1 + 0{VlnV) (59) 

Note that Eq. ( |59| ) only holds at intermediate distances from the origine, since, on the one 
hand, n and m have to be large enough in order to apply th stationary phase method, but, on 
the other hand, the argument of Kq in Eq. ( p8| ) should be small. However, as shown in Ref . |3Q| , 
the asymptotical result in Eq. (p9|) serves as a very good estimate when n and m are sufficiently 
large (in practical, as soon as {m? + m^)^/^ > 3). Since we are also interested to know /i±i,o and 
^o,±i precisely, we list below some particular values of P{n, m\C,) in the immediate vicinity of the 
tracer. Exact values of P{n, m\C, — *■ 1) are presented in Ref.[Q. From these, we find 



V„Fo,o = -1 + 0{Vln{V)), VnFifl = -3 + S/tt + 0{Vln{V)), A„Fi,o = 2 - S/tt + 0{Vln{r)) 

(60) 

V„i^_2,o = 3 - 8/7r + OiVlniV)), V„i^_i,o = 1 + 0{Vln{V)), A^F.^^o = 2 - S/tt + 0{Vln{V)) 

(61) 

V„F_i,i = -1 + A/t: + 0{Vln{V)), V„i^o,i = 1 - 4/7r + 0{Vln(V)), A,„Foa - -2 + 8/n + 0(Vln{V)), 

(62) 

which allow us to determine the local deviations in the immediate vicinity of the tracer. We find 
that in the limit (c) the local deviations are to the leading order in V: 

hi,o ^ apLp + 0{P^) h^,,o = -apLp + 0{P^) ho,i^O{P^), (63) 

where a is a numerical constant, 

P(2,0|l)-P(0,0|l) 7T-2_ 

"^' 4 + P(2,0|l)-P(0,0|l) ^^""-^'^ ('^) 



Finally, inserting Eqs.(63) to Eq.(||), we find the following expression for the tracer terminal 
velocity in the limit (c): 



_ il-pL)Do 

Vtr — 



l-^+0(i^f)\ (65) 



acrpL 

Note, that similarly to Eq.(|5l|), in the limit (c) the tracer velocity Vtr vanishes in proportion 
to Dq; the prefactor is, however, different from that in Eqs.(Ma) and Eq. (pi]) . The Peclet number 
in this case is again P ~ (1 — pl)/pl, which means, in view of the consistency with the initial 
assumption, that the limit (c) is only realized and Eq.(|65|) is only vaHd when p^ ~ 1. We also 
remark, that a similar to Eq.(pq) expression for the tracer velocity with, however, a bit different 
value of a, has been obtained prior in pG], which studied the stationary velocity of a driven tracer 
particle in a two-dimensional lattice gas with a conserved particles number. 
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Lastly, using Eq. (BO) and Eq. (p9[), we find that at a finite, but large enough distance from 
the tracer, 

(l-p)(l + a) , ( {n + lf + m' \ 

which is again a non-separable function of variables n and to. Profiles defined by Eq. ( |66[) are 
depicted in Fig. 2c. 

5.2 Limit of a large particles diffusivity. 

We turn next to the analysis of the asymptotical behavior of Vtr in the limit of very fast monolayer 
particles diffusion supposing that both / and g are non-zero, which will be called the limit (d). 
Clearly, when Dq = oo, one should recover for Vtr the trivial mean- field result in Eq.(|^). We 
search, hence, for the first correction to this result, representing the tracer velocity as 

Now, one can readily notice that in this case we deal, in essence, with a situation which is quite 
similar to the limit (c) studied in the previous subsection. Namely, both S and V are small and 
hence, ( should be of order of unity, ( ^ 1. Consequently, all the analysis of the last subsection 
applies here, except that C takes a bit different form. 

Noticing next that here ( ^ 1 — S and calculating the integral in Eq. ( p2[ ) along essentially the 
same lines as it was done in the previous subsection, we find that in the limit of fast particles 
diffusion P(n,jn\(^ —^ 1) is given by : 

P(n,m|C -> 1) ln[^) - -^ ^ - ^ + 0{^-^), (68) 

which expansion holds for large enough values of the variables n and to. Consequently, we find 
that in the limit (d) the local deviations from the equilibrium mean density in the immediate 
vicinity of the tracer obey Eqs.(|63[) and hence, the first correction to Eq.(||) reads 

OL(7 

^(i)--Pl(I-Pl)— (69) 

As it could be expected intuitively, V(i) appears to be negative, i.e. the actual tracer terminal 
velocity is lower than that given by Eq.(^). Note also that V(i) is a non-monotonous function of 



PL and is maximal when p^ = 1/2. We finally remind that the expansions in Eq.(67) and Eq.(p£ 
hold, by definition, in the limit Do —f cx3, which means they only apply when both V and 6 are 
much less than unity. 

We remark that Eq.(|69|) and Eqs.(|6^) still hold when /, g — 0. On the other hand, it will be 
shown in the Appendix that in the lim,it (d) the amplitudes and the characteristic lengths of the 
density profiles attain different forms depending whether 6 is equal to zero or not. 

Lastly, the form of the density profiles at a finite distance from the tracer is very similar to the 
one found in the limit (c). We have that, at finite but large enough distance, ft.„.m is given by: 

. - P(l-p)(l + a).^^ /(n+2)!±4\ (70) 



47r tDq \{n- 1)^ -|- to' 

The remarks following Eq.(|6q) still hold. 
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6 Asymptotic forms of the density profiles at large dis- 
tances from the tracer. 

In this section we consider the asymptotical behavior of the density profiles at large distances 
from the stationary moving tracer particle. The first subsection deals with the simple special 
case Do — 0. The second subsection will be devoted to the analysis of the asymptotic forms of 
hnfi in the limit n — > cx), which will be based on our Eqs.(pO[) to (p2). Further on, in the next 
subsection we will examine a more complicated question of the decay of hn,o at large separations 
past the tracer, where we again discuss separately the behavior in the cases with non-conserved and 
conserved particles numbers. Actually, we observe here a very spectacular effect of a qualitative 
change of the asymptotical behavior; the density relaxation to the average value pL proceeds in 
the former case exponentially, while in the latter case it is described by a slow algebraic function of 
the distance. Lastly, we will study the forms of some integral characteristic of the density profile, 
such as, for instance, the global deviation of the density from the equilibrium value in the domains 
n > 1 and n < — 1. 

Asymptotic behavior of /i„_o as n — > ±cxd can be most conveniently studied if we introduce the 
generating function of /i„^o of the form 

/i(z)=^Vo^" (71) 

— oo 

Setting then ttt, — in Eq.(pO|) and summing over all /?,, we find that h{z) is given explicitly by 



z{z~l)({l+V)h^,o + VpL)+{l~z){h^i.o-VpL) liz-z,)iz-z,) ^_^ 
hyz) = , "oia/t 77 r (72) 

^{z^Zi){z-Z2){z-Z3){z-Z,) ' \l {Z ^ Z,){Z - Zi) 

where the roots z, are as follows: 



zi = 



3. = -'--:- ' -,; r"-^-M '-(i.n 



Zi 



i + Vi + T' + C 


c 


1 + Vi + V-C 


c 


l + Vl + T'-C 


c 


i + Vi + ^ + C 



'/l + VT+V + C 



y-(l + V); 



.3 = c + \/i 1 ) -(1+n 



■1 + VT+^ + C 

c 



y-{l+V) (73) 



and obey 



< zi < Z2 < 1 < Z3 < Z4, (74) 



such that h(z) is analytic in the annular region of inner radius Z2 and outer radius z^. This implies 
that both zi and Z4 are irrelevant to the large scale behavior and the dominant contribution to 
the asymptotic form of hn,o &s n —^ ±(X) comes from the the behavior of h{z) as z — > Z2 or z — > Z3. 
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More precisely, the asymptotic form /i„^o as n —s- +00 stems from the behavior of h{z) in the 
vicinity of z = Z3, while hnfl as n — > — cx) is associated with the corresponding behavior near 

Z = Z2- 

6.1 Density profiles in the special case Dq = 0. 

First of all, we study the special case Dq = 0, for which the roots zi and Z2 collapse to the same 
value (1 — Pl)/{^ + if + g)T/T*). In this special case, which corresponds to the original Langmuir 
model, it is possible to determine the density profile exactly. Indeed, here Eq. (f2|) reduces to 

h[z) = hi^o — 



z(l - Z2/2;) 

^ CSO 



(1 - Z2)Z. 



n=0 



,,,(. + 2 li^iilf^), ,75) 



which implies 



where h\x) and /i_i.o are given to the leading order by Eqs.(|70[) and ( |44| ) respectively. These 
expressions confirm the step shaped profile in front of the tracer particle (as anticipated in section 



.„,o = 0, for n > 1 , /i„,o = /i-i,o ( ^ + ^^ + ^^^1^ ) , forn < 0, (76) 



5.1), such that the only perturbed sites are the sites visited by the tracer particle, and the site just 
in front of it. These expressions can be also readily recovered by considering the limiting form of 
the equations (|o|)- (^) when Z^o = 0. 

6.2 Asymptotic forms of the density profiles at large distances in front 
of the stationary moving tracer. 

From now on, we assume that the diffusion coefficient Dq is not equal to zero, and we deduce 
the asymptotical behavior of /i„^o in the limit n — > ±00 from the analysis of the singularities of 
the generating function h[z). One has that in the vicinity of z = Z3 the generating function /i(z) 
obeys 



Z3(Z3 - 1) ((1 + V)h^,a + ^Pi) + (1 - 23) (/l-l,0 - T'pi) 



KA = ^ ,, ^, ^, ^^^+0(v/i^^) (77) 

\/{z-i- ZY)\Z-i- Z2)\Zi- Z-i) \JZ-},-Z 



Next, following Darboux method |27| ] or singularity analysis of generating functions |31|, we find 
that /i„^o in the limit n — > 00 follows 

Z3(z3-1)((1+7')/ii,o + 7'pl)+(1-Z3)(/i-i,o-7'pl) \ 1 
^Z3 (z3 - Z2) (z4 - Z3) (z3 - zi) ^J^mz■^ 

For notational convenience, the latter expression can be rewritten as 
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hn,0 '^ K+- 
where the decay amplitude is given by 



exp I — n/A+ j 

771/2 



(79) 



(^3 - 1) ((1 + 7')/7i,o + Vpl) + (1/^3 - 1) (/i-i.o - Vpl) 



K.^ 



while the characteristic decay length obeys 

A, 



/i + Vi + v-cy 


-{l + V) 




(80) 



-1/4 



l/ln(z3) 



(81) 



Therefore, the density in front of the stationary moving tracer approaches the equilibrium value 
PL exponentially with the distance. 

Note now that the decay amplitude K+ and the characteristic length A+, Eqs.(|8^) and (|8l|), 
depend on C and on /i±i,o- It may be thus instructive to analyze the asymptotical forms of K^ and 
A+, using the explicit results for /i±i,o obtained in the previous section. This analysis is presented 
in Appendix 2. 

6.3 Asymptotic forms of the density profiles at large distances past the 
stationary moving tracer. 

We first note that one of the roots of the generating function, namely Z2 , gets equal to unity when 
both / and g are strictly equal to zero, which results in the exact cancellation of the multiplier 
1 — 2:3 both in the nominator and the denominator in Eq.(^I]) . This shows that in the limit 
when exchanges with the reservoir are forbidden, qualitative changes in the singular behavior of 
the generating function at the vicinity of 22 will appear as compared to the case when f,g>0. 
Consequently, we have to consider separately the behavior in the case of non-conserved particles 
number, when exchanges with the reservoir persist, and the case when both / and g are equal to 
zero while their ratio is kept fixed. 

A. Non- conserved particles number. From Eq. (|7l|) one finds that in the vicinity of z = Z2 the 
generating function h(z) behaves as 

y/ (Z2 - Zl) (Z3 - Z2) (Z4 - Z2) ^/Z-Z2 

which implies the following form of /7„_o in the limit n — *■ — c». 



exp I 



/i„,o ^ K^ - 



with 



K^ 



(-n)-i/2 ' 
(22 -!)((!+ ^)/^i,o + ^Pl) + (1/^2 - 1) (/^-i,o - ^Pl) 



(83) 



87r 



Vl + v 


C 


-^)'- 


-(l + T') 
(84) 



■1/4 
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and 

A_ = -l/ln{z2) (85) 

B. Conserved particles number. Suppose now that both / and g are equal to zero, i.e. S — 
0, while their ratio is fixed and given by f/g = Pl/(1 ^ Pl)- As we have already remarked, 
this situation corresponds to the usual model of a two-dimensional hard-core lattice gas without 
exchanges with a reservoir. Here, the generating function in the vicinity of its singular point 
z — Z2 ^ I obeys 

1 



Hz) 



2PpL + (1+^)^1.0-^-1.0 



^0,l\/^3 - 1 



^^=VI^ + 0((z-l)V2), 

(86) 



which yields the following result for /i„.o in the limit n 



where 

Remarkably enough, in this case the correlations between the tracer position and the particles 
distribution vanish only algebraically slow with the distance! This implies, in turn, that in the 
conserved particles number case, the mixing of the monolayer is not efficient and the medium 
"remembers" the passage of the tracer for a long time which signifies strong memory effects. 
We note also that the algebraic decay of correlations in this model has been predicted earlier 
in Ref.|2^]. However, the decay exponent has been erroneously predicted to be equal to 1/2, as 
opposed to the value 3/2 given by Eq. (|8^) . As well, the amplitude K- happens to have a different 
sign, compared to that in Ref. [pOJ , which invalidates the conclusion that the overall relaxation to 
the equilibrium value p^ is a non-monotonic function of the distance. 

Asymptotic forms of K^ and A_ in different limiting cases are presented in Appendix 2. 

6.4 Integral characteristics of the density profiles 

First of all, we address the question whether the driven tracer, which induces an inhomogeneous 
density distribution in the monolayer, shifts the equilibrium between adsorption and desorption, 
i.e. whether it changes effectively the equilibrium density in the monolayer. The answer is trivially 
"no" in the case when the particles number is explicitly conserved, but in the general case with 
arbitrary / and g this is not at all evident: desorption events are certainly favored in front of the 
tracer, while the adsorption events are evidently suppressed by the excess density. On the other 
hand, past the tracer desorption is diminished due to the particles depletion while adsorption may 
proceed more readily due to the same reason. It is thus not at all clear a priori whether these 
two effects can compensate each other exactly, in view of the asymmetry of the density profiles. 

For this purpose, we study first the behavior of the integral deviation 51 of the density from 
the equilibrium value p^, i.e. 

n— — oo m— — cxD 
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As a matter of fact, this integral characteristic can be computed straightforwardly from Eq.(26). 
Setting z and w equal to unity, we find that even in the case when (5 7^ 0, the integral deviation 
n = 0. This means that the inhomogeneity of the density distribution in the monolayer created 
by the tracer does not perturb the global balance between the adsorption and desorption events. 



An analogous result has been obtained for the one dimensional problem in Ref.|21|. 

Next we would like to check whether such a compensation holds on the axis of the tracer 
motion, i.e. whether the integral deviation ^x along the X-axis is zero or not. Setting z = 1 in 
Eq. ( [Til ) , we find then the following general result 



n— — 00 



^x= Y. ''"^o = "\/tt^ ''•'''' ^^°^ 



which implies that such an exact compensation is not realized on the X-axis only, except for the 
case i5 = 0, which also does not seem to be a trivial fact. Moreover, it follows from Eq.(|9y) that 
f2x is always negative, provided that ft,o,i is positive definite, which means apparently that on 
the X-axis the equilibrium is shifted towards desorption. Note, however, that fix always remains 
small and tends to zero in the limits Dq — j and Dq — > cx); hence, f2jf should be a bell-shaped 
function of the particles diffusivity. 

Finally, we look at the total excess "weight" of the condensed region in front of the tracer 
particle, defined as 

00 
fix + = Y^ /i«,o (91) 

Using the expression for /i„^m in Eq. (pO[) , we find for r2x+ the following asymptotical results: 

(a) In the limit of small particles diffusivity and (5 7^ (T' ::^ 1), 

(-, Pl(X- PL) , . 

This quantity is smaller than one, but finite in the general case. Moreover, it appears to be non- 
monotonic with respect to p^. Lastly, it tends to zero if r* tends to zero, i.e. when the exchanges 
with the vapour phase proceed infinitely fast. 

(b) In the limit of small particles diffusivity, (5 = and "P <C 1, f2x+ ^ 1, which result can be 



obtained from the previous one (92) by taking the limit (5^0, and recollecting that p^ <C 1 in 

the limit (b). 

(c)In the limit of small particles diffusivity, (5 = and V ^\., 

(l + a)(l-pi) 



(t^)- '») 



f2x-i- ~ In 

an M — PL 

Since here pL ^ 1, i^x+ appears to be small. 

(d) In the limit of fast particles diffusion {V 3> 1), i^x+ approaches a constant value 

(1 + a)a^ PLJl ~ Pl) , ( DqT \ 

^x+ 7^ '■^[^1 T • (^4) 

ttDot V(t^(1 -pLJ/ 

Note that ^x+ is non-monotonic with respect to p^, which generalizes the result announced in 
pO| to the case with non-conserved particles number. Note also that ^x+ ^ as Dq -^ 00, which 
shows that the condensed region disappears if the diffusion processes become efficient enough to 
mix the adsorbed monolayer. 
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7 Conclusions 

To conclude, we have studied analytically dynamics of a driven probe molecule in a two-dimensional 
adsorbed monolayer composed of mobile, hard-core particles undergoing continuous exchanges 
with the vapour phase. Our analytical approach was based on the master equation, describing 
the time evolution of the system, which allowed us to evaluate a system of coupled dynamical 
equations for the tracer particle velocity and a hierarchy of correlation functions. To solve these 
coupled equations, we have invoked an approximate closure scheme based on the decomposition of 
the third-order correlation functions into a product of pairwise correlations, which has been first 
introduced in Ref . |23| ] for a related model of driven tracer dynamics in a one-dimensional lattice 
gas with conserved particles number. Within the framework of this approximation, we have 
derived a system of coupled, discrete-space equations describing evolution of the density profiles, 
as seen from the moving probe, and its velocity Vtr- We have shown then that Vtr depends on the 
monolayer particles density in front of the tracer, which is itself dependent on the magnitude of 
the velocity, as well as on the rate of the adsorption/dcsorption processes and the rate at which 
the particles can diffuse away of the tracer. As a consequence of such a non-linear coupling, in 
the general case, (i.e. for arbitrary adsorption/desorption rates and particles diffusion coefficient), 
Vtr has been found only implicitly, as the solution of a certain non-linear equation relating its 
value to the system parameters. This equation simplifies considerably in the limit of small or large 
particles diffusivity, in which two cases the tracer velocity is calculated explicitly. Further on, we 
have found that the density profile around the tracer becomes strongly inhomogeneous: the local 
density of the monolayer particles in front of the tracer is higher than the average and approaches 
the average value as an exponential function of the distance from the tracer. The characteristic 
length and the amplitude of the density relaxation function are calculated explicitly. On the other 
hand, past the tracer the local density is lower than the average; we show that depending on the 
condition whether the number of particles in the monolayer is explicitly conserved or not, the local 
density past the tracer may tend to the average value either as an exponential or as an algebraic 
function of the distance, revealing in the latter case especially strong memory effects and strong 
correlations between the particle distribution in the monolayer and the tracer position. 
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8 Appendix 1: Series representation of F^^? Eq. (pi| ), in 
the limit (a) 

We calculate to the first order in powers of Dqt* /a^ the functions Fn.m involved in Eq.(|39|): 



Fo,o^l + 2— ^^+o((^_)^), (95) 

r*fy(o)+a(/ + 5)/T*J '^ 
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T*fy(0)+a(/ + .g)/r*j '^ \ a' J 



'(v(")+a{f + g)/T*^ '^ 

^-1.0-7 \- ^ ; Ts 2~ + *^(( — 2~) ' 

{vW+a{f + g)/r*) t*^{v(o) + aif + g)/T*) " ^ ^ 

(98) 
and 

y(0)2 
-F-2.0 



{v(")+aif + g)/T*) 

aT/(") ((/ + ,g)2(T/(i)a/T* + {a/r*)^) + (/ + 5)F(")(F(i) - 2a/T*) - F^^) ^ * 
-I- 2 -^ - - I- 

r*(y(")+a(/ + .g)/T*j '^ 

+ 0((^)^), (99) 



while 

i^2,o = i^i,±i = i^o,2 = o((:^)2) (100) 

9 Appendix 2: Explicit results for the amplitudes and the 
decay lengths of the density profiles at large distances 

9.1 In front of the tracer 

We find the following explicit asymptotical forms for the decay length and the amplitude: 
(a) For smaU Dq, 5 > and T' > 1, 

+ V DoT*{p+{f + g)T/T*) J' ^ ^ 

and 

Pl{1 - PL){.f + gf/' (1 + (/ + 5)t/t*)3/2 . a2 x3/2 



K, 



20F / /. N / A^/HDot 



{pL + i.f + 9)T/T*y 



i^y ■ <'-> 



which means that A+ is logarithmically small with Do, while K+ is large. Note also that K+ can 
be a non-monotonous function of p^ when (/ + g)T Jt* ^ pj^. 
(b) For smaU Dq, 5 = and P > 1, 

^-s^-"-'(^)' <-' 
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and 



-p3/2 
K+ ^ — =-, (104) 

4v27r 



i.e., A+ is logarithmically small with V, while ii'+ is large. 
(c) For small Dq, S ^ and T' < 1, 



A+ ^ 1/V = -^^^; (105) 

1 -PL 



and 



X+^l(l + a)p,7'^/2 = i(l + a)il-^ (106) 

Hence, in this case A4. is large, since it is inversely proportional to the small parameter V, but 

K+ is small. 

(d) For Do —> 00 and 6 > 0, 



' -^-\7?^ (107) 



^ 2V6 \{f + 9) 

and 

K -.n , .,^„_-n/x/„2^l/^_ (l + «K(l-PL)aV'^'(/ + 5)^l/4 



'-+ 



a.oK.(v.')-^ ^'-°rX"" ("i^) <-' 



i.e., here, likewise to the case (c), A+ appears to be large and K^ is small. On the other hand, if 
we suppose that S — (no adsorption/desorption) and Dq — > 00, we find a bit different forms of 
K^ and A+; namely, 

A+ - :^ = jr^^ (109) 

and 

j^^_ (l + a)pLV'/^ _ (1 + c.)pl / ^^(1-Pl) \3/2^ ^^^^^ 

i.e., similarly to the above considered case, A+ is large and K^- is small. Note also that, despite 
the fact that the tracer velocity in the cases 6 > and (5 = is given by the same expression in 
Eg. (|69[) , the characteristic properties of the density profile appear to be different when S is zero 
or S ^0. 

9.2 Past the tracer 

A. Non-conserved particles number. 

Consider next the asymptotical behavior of K- and A_ in the limit of slow and fast particles 
diffusion. 
For small Do, S > and V ^ I, which corresponds to the limit (a) in our previous notations, 



VI + if + q)t T* J 
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and 

^ PL ^ (/ + g)(l + (/ + g)r/r*) ^ a^ y/2 

i.e. A_ tends to a constant value when I?o — > 0, which means that here A_ > A+, while \K^\ (x 

1/9 

Dq ' diverges when Dq -^ 0. Note, however, that we evidently have \K-\ <C K+. This means 
that in the limit (a) the density profiles around the tracer are strongly asymmetric; in front of the 
tracer we have a condensed, "traffic" -jam-like region, characterized by a high amplitude but of a 
relatively short spatial extent, while past the tracer there is a depleted region which extends on 
much longer scales but has a considerably smaller amplitude. 

Next, for high particles diffusivity and (5 > 0, which corresponds to the limit (d), one encounters 
an exactly opposite situation. Here the values of the characteristic lengths and the amplitudes of 
the condensed and the depleted regions almost coincide, i.e. K- w — i4r+ and A- w A+, where K+ 



and A-|_ are given by Eqs. (|lO^ ) and (107). 

B. Conserved particles number. 

Consider next the asymptotic forms of the decay amplitudes in Eq. (|8g) . Using the results of 
the previous section, we find then the following asymptotical results: 
In the limit of small Dq, S = and V ^ 1, which corresponds to the limit (b) of Section 5, 

which means that the amplitude is large in this limit, since p/^ <C 1. 
In the limit of small Dq, (5 = and T' <C 1, (limit (c)), 



Since this limit can be only realized at sufficiently high particles densities, pL ^ 1, we have then 
that in this case the amplitude should be small. 
Lastly, in the limit Da —t oo (limit (d)) we find 

2 V ttDot / ' ^ ' 

i.e. which signifies that in this case the decay amplitude is small. 
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Fig.l. Two-dimensional lattice of adsorption sites partially occupied by identical, mobile 

hard-core particles (grey spheres) undergoing continuous exchanges with the reservoir - the 

vapor phase. Particles desorption and adsorption probabilities are denoted by g and /, 

respectively. The dark grey sphere with an arrow denotes the tracer particle, whose motion is 

completely directed by external field E. 
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Fig. 2a. Sketch of the density profiles at finite distances from the tracer. The black circle denotes 

the tracer particle and the arrow represents the direction of the external force exerted on the 

tracer. Limit (a): Dq -^ 0,6 ^ 1, and V ^ 1. 
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Fig. 2b. Sketch of the density profiles at finite distances from the tracer. Limit (b): 

Do->0,S = 0, and-p > f. 
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Fig. 2c. Sketch of the density profiles at finite distances from the tracer Limit (c): Dq ^ 0, (5 = 0, 

and P < f . 
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